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INTRODUCTION 


The use of computational methods for the prediction of multi-dimensional recir- 
culating flows has been continuously increasing over the recent years. As the nu- 
merical techniques become more powerful, they are being applied to even more chal- 
lenging problems encountered in combustion chambers, gas turbines, rotating machin- 
ery, heat exchangers, and other devices. Although computation is far less expensive 
than full-scale testing, the cost of a computational run is still substantial. 
Therefore, attempts are continually being made for improving the accuracy and effi- 
ciency of numerical techniques so that the predictions of a given accuracy can be 
obtained at a modest cost. 

The central feature in any flow prediction method is the treatment of the cou- 
pling between the momentum and continuity equations. In natural-convection flows, 
the energy equation also becomes strongly coupled with the momentum equations. Be- 
cause of the nonlinear nature of the coupling, these equations are solved iterative- 
ly. Iterative methods are often prone to slow convergence, divergence, and extreme 
sensitivity to underrelaxation factors. 

The aim of the present research is to develop more efficient and reliable solu- 
tion schemes for the coupled flow equations. Such schemes will significantly reduce 
the expense of computing complex flows encountered in combustion chambers, gas tur- 
bines, heat exchangers, and other practical equipment. 

In the work completed so far, a technique employing norm reduction in conjunc- 
tine with the successive-substitution and Newton-Raphson techniques has been devel- 
oped. Also, a block-correction procedure for the flow equations is currently being 
formulated and tested. 


NORM REDUCTION TECHNIQUES 

The development of a number of methods for solving strongly coupled equations 
has been reported in reference 1. The recommended method there is a combination of 
the successive-substitution and Newton-Raphson methods coupled with a norm reduction 
technique. 

The central idea of the method will now be described. The linearized discreti- 
zation equations are solved by a direct method, such as the sparse-matrix LU decom- 
position. The linearization can be of two kinds. If the unknown coefficients are 
simply evaluated from the currently available values of the dependent variables, the 
linearization is called successive substitution. In the Newton-Raphson method, the 
anticipated change in the coefficients is taken into account via their first deriva- 
tives with respect to the dependent variables. 


* Work done under NASA Grant NAG 3-596 
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The solution of the linearized equations predicts a change in the values of the 
variables. However, instead of accepting this change as is, it is multiplied by a 
constant. (This constant can be imagined to be a kind of under- or over-relaxation 
factor.) The value of this constant multiplier is found by requiring that the norm 
of the residual vector (i.e., the sum of the squares of the residuals of all the 
equations) be a minimum. This minimization search produces a kind of "optimum" 
underrelaxation of the dependent variables. 

The norm minimization (or reduction) can be used with either successive-substi- 
tution or Newton-Raphson linearization. When the current estimates are close to the 
final solution, the Newton-Raphson technique is very efficient. But for initial 
guesses that are far from the solution, the technique often diverges. Therefore, 
the practice found satisfactory for all the flows tested was to employ the succes- 
sive-substitution linearization until the norm became less than a small quantity and 
then to switch to the Newton-Raphson linearization until the final convergence. Be- 
cause of the combination of the two linearization practices, the technique is called 
the hybrid method in reference 1. 

The hybrid method and many other alternative schemes (some of which are based 
on the Broyden methods described in reference 2) were tested on two flow configura- 
tions: the flow in a driven cavity at different Reynolds numbers and the natural 

convection flow in an enclosure with hot and cold walls at different Rayleigh num- 
bers. 

For the driven cavity problem, solutions were obtained by the hybrid method in 
at least one-third the computer time required for the iterative method SIMPLER. For 
the natural convection problem, most methods either diverged or converged extremely 
slowly as the Rayleigh number was increased. The hybrid method, however, converged 
rapidly and required a very modest amount of computer time. By using the hybrid 
method, it was possible to obtain convergence in 25 iterations at a Rayleigh number 
of lO? when a zero initial guess was used for all variables. If the results for a 
lower Rayleigh number could be used as the initial guess, it was possible to obtain 
solutions at a Rayleigh number of 10^. It is believed that, for the first time, 
solutions have been obtained for such a high Rayleigh number (with a Prandtl number 
of 0.71). Iterative methods such as SIMPLER failed to converge even after 1000 
iterations . 


BLOCK CORRECTION TECHNIQUE 

Work is currently in progress on another approach for accelerating the conver- 
gence rate of an iterative procedure such as SIMPLER. In the proposed technique, 
the velocity and pressure values are adjusted through a block correction procedure. 
Here the calculation domain is considered to be composed of several large blocks. 

Each block contains a number of grid points (and hence control volumes). It is pro- 
posed that the values of a variable for grid points within a given block will re- 
ceive a uniform correction. These corrections are calculated such that the integral 
conservation of momentum and continuity is satisfied for each block. The block- 
correction. equations thus resemble the discretization equations for momentum and 
continuity but are formulated on much coarser grid. The solution of these equations 
by a direct or iterative method is rather straightforward. 

Initial testing of this apprach shows that, for fine grids, there is a notice- 
able improvement in the computational effort required to obtain a converged solution. 
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